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V 


At  the  initial  stage  of  development  of  wind  waves  in  the  sea,  the  energy  dis- 
sipation is  determined  by  the  molecular  viscosity.  In  this  case,  the  dissipated 
energy  averaged  over  the  wave  period  may  be  calculated  from  the  formula  (see  for 
example  Ref.  1):  c 6 f. 


/773 


Emolec  ~ 2pvk’cla1, 


(1) 


where  p is  the  density  of  the  liquid,  V is  the  kinematic  coefficient  of  molecular 
2tt 


viscosity,  k = 


X is  the  wavelength,  c is  the  wave  propagation  velocity,  and  a 


is  the  wave  amplitude.  As  the  wind  waves  develop  further,  the  principal  role  in 
their  energy  change  due  to  scattering  begins  to  be  played  by  turbulence.  The  mag- 
nitude of  the  energy  dissipated  in  turbulent  wave  motion  may  be  calculated  from 
the  formula 


where 

V 

(2) 

u ---■  ao>*‘  ~kt  cos  {kx  ■ ~ oj t). 

(3) 

The  bar  in  Eq . (2)  denotes  averaging  over  the  wave  period,  and  V is  the  volume. 
For  the  coefficient  of  turbulent  viscosity  vturj}»  t^ie  following  expression  common 

in* semiempirical  theory  may  be  used  in  Eq.  (2): 


du 

Vrt  ==  ’ ~di. 


(*) 


where  SL  is  a quantity  representing  the  turbulence  scale,  or  the  Prnndtl  mixing  length. 

It  is  convenient  for  later  purposes  to  introduce  the  coefficient  of  integral 
turbulent  viscosity  vturjj»  given  by  the  relation 

Emolec  ~ ,'7„r()  *"or  v 3 v i„<b  • ^ 

Tlius,  to  calculate  the  dissipation  energy,  it  is  necessary  to  know  the  expres- 


sion for  the  turbulence  coefficient  V 


turb 


Since  turbulent  viscosity  characterizes 


the  properties  of  turbulent  motion  in  a wave,  v ^ should  be  determined  by  all  the 

wave  parameters  (wave  height  h = 2a,  wavelength  X and  wave  period  T) . It  is  easy 
to  see  that  these  characteristics  are  exhaustive,  since  the  integral  turbulent  vis- 
cosity V .is  independent,  of  the  coordinates,  and  one  can  therefore  write 
turb 


vtufb  -Hh.X.T). 


(6) 


Numbers  in  the  right  margin  indicate  pagination  in  the  original  text. 
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r*  -w  «•» 
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Considering  the  dimensions  of  the  quantities  in  Eq.  (6)  and  denoting  by  6 the 
dimensionless  parameter  ^ representing  the  wave  steepness,  we  can  obtain  the  most 
general  relation  for  the  coefficient  vtur^: 

Concerning  the  type  of  the  function  'l'(fi)  in  Eq.  (7),  the  following  statement 
may  be  made:  the  steeper  the  waves  (the  larger  6),  the  greater  the  velocity  gradients 

and  energy  dissipation.  Therefore  \>turb  and  hence,  i^(6)  should  increase  with  in- 
creasing 6,  i.e. , 

(jr* / V W>°-  j (8) 

^ T M's-:'1 

In  order  to  be  able  to  discuss* various  theoretical  models  of  turbulence  in 
wave  motion  and  compare* them  with  existing  indirect  and  empirical  dataKwe  will 
postulate  that  the  functional  relation  (6)  may  be  written  in  the  form'af  a power 
law:  i 

Vt„b  "h’XYrP*  (9> 


where  n is  a dimensionless  constant  coefficient.  On  the  basis  of  dimensional  con- 
siderations, it  follows  from  (9)  that 

a -|-  v --  2;  p — 1.  (10) 


Using  the  well-known  relations  of  wave  theory 


T ” C; 


and  introducing  the  notation 

P 1=3  w ’ 


,-37  ’ 


K =-•  2nn, 


where  W is  the  velocity  of  the  wind  which  generated  the  waves,  considering  Eqs.  10, 
we  can  write  Eq.  (9)  in  the  form 


k^T. 


If  use  is  made  of  the  experimental  fact  that 

8 — / (P). 

one  can  represent  Eq.  (13)  in  the  form 

W* 

W -K-f*®' 


where  <j)(3)  will  depend  on  both  the  type  of  relation  (14)  and  on  the  magnitude  of  the 
exponent  a in  Eq.  (13).  Relation  (15)  obtained  shows  that  the  turbulence  coefficient 
depends  on  the  wind  velocity  W and  on  the  stage  of  wave  development  3 (3  is  frequently 
called  the  wave  age).  The  latter  fact  leads  to  the  important  conclusion  that  the 
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question  of  selection  of  the  exponent  a is  not  a formal  one,  since  a represents 
the  degree  of  influence  of  the  scale  of  wind-driven  waves  on  turbulent  exchange  in 
the  upper  layer  of  the  sea.  It  follows  from  Eq.  (15)  that  in  the  case  of  fully 
developed  waves  (8  - 1),  the  turbulence  coefficient  depends  only  on  the  wind  velo- 
city. In  contrast  to  the  previously  known  Ekman-Rossby2  quadratic  dependence 

Vturb^’  We  °bta*ned  a cubic  one.  In  Ref.  3,  such  a dependence  was  obtained 

empirically.  The  data  given  in  Ref.  4 may  also  be  regarded  as  confirmation  of  the 

cubic  dependence  vturjj(w)* 

Since,  as  has  been  shown  experimentally,  the  wave  height  and  length  increase 
with  increasing  8,  V should  also  increase  with  8,  i.e., 

(?0>vx„w  (?i)  for  P»>P»-  (16) 


If  the  form  of  function  (14)  is  considered  to  be  reliably  established  empirically, 
condition  (16)  may  serve  as  a criterion  for  selecting  a.  On  the  other  hand,  as  we 
have  noted  above,  the  selection  of  a in  Eq.  (13)  should  also  be  consistent  with  the 
most  general  dependences  (7),  (8). 

We  will  go  directly  to  the  determination  of  the  function  4>(3)  in  Eq.  (15)  for 
different  values  of  a in  Eq.  (13).  We  will  use  the  relationship,  known  from  field 
observations,  between  the  wave  steepness  6 and  wave  age  8»  Such  a relationship 
was  first  studied  Sverdrup  and  Munk.5  This  empirical  relationship  was  later  refined 
by  Neumann^  and  may  be  written  as  follows  (Neumann's  data,  obtained  for  so-called 
"substantial"  waves,  were  converted  for  medium  waves  according  to  statistical 
wave  theory®) : 

1 

S -0,062  for?  < '3  . 

8 =-  0,1073«— l,W7®  for  y < P < *•  (17) 


Recently,  Krylov^  obtained  a theoretical  dependence  of  the  wave  steepness  6 on 
8,  written  in  the  same  form  (we  use  the  rounded  values  given  in  Ref.  7):  /77.‘j 

8 - 0,06  for  ? sS  0,375, 

0.0225  _ „ , (IS) 

8 — --*• — for  0,375<P<1. 

Figure  1 shows  both  of  these  curves  along  with  the  experimental  data  of  Sverdrup 
and  Munk,  Neumann,  and  Vilenskiy  and  C.lukhovskiy It  is  evident  that  the  agreement 
of  the  curves  is  completely  satisfactory,  and  therefore  either  one  can  be  used  for 
our  purposes.  We  shall  subsequently  use  relation  (18),  since  it  was  obtained  theore- 
tically and  is  written  in  the  form  most  convenient  for  calculations.  In  Fig.  1,  the 

region  least  studied  experimentally  Is  the  range  of  8 values  from  0 to  0.4  (more 

accurately,  from  8 , = — — , where  c , = 0.23  m/scc;  calculations  show  that  8 , 

min  W min  min 

is  close  to  zero).  Some  experimental  results  indicate  that  when  8 < 0.4,  the  wave 
steepness  increases  sharply.  To  simplify  further  calculations,  we  will  con- 
sider the  waves  in  the  interval  0 < 8 < 0.4  to  be  turbulent  ones  with  a constant 
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Fig.  1.  Relationship  between  wave  steepness  Fig.  2.  Turbulent  viscosity  coefficient 
6 and  wave  age  3:  ^ vs  wave  age  3: 

1 - Neumann’s  empirical  curve,  ^ 1 - curves  corresponding  to  different 

2 - Krylov's  theoretical  curve;  values  of  a in  Eq.  (13),  relation 

3 - empirical  data  of  Vilenskiy  and  (19)  being  used; 

Glukhovskiy ;®  2 - Neumann's  empirical  curve. ^ 

4 - empirical  data  of  Sverdrup  and 

Munk.5 


steepness.  Since  the  case  of  developed  waves  (3  > 0.4)  is  of  principal  interest 
in  determining  the  turbulent  viscosity  in  the  upper  layer  of  the  sea,  this  assump- 
tion will  not  appreciably  affect  the  results  obtained.  In  addition  to  relations 
(17)  and  (18),  there  are  also  the  data  of  Darbyshire^jlO  pertaining  to  the  dependence 
of  the  wave  steepness  6 on  3.  The  structure  of  the  relationship  6(3)  in  them  is 
the  same  as  in  relation  (18);  therefore,  in  determining  the  function  vturj,^’ 

use  of  Derbyshire's  data  will  give  the  same  results  as  when  formula  (18)  is  employed. 
The  use  of  the  universal  relation  6 = f(3)  instead  of  all  the  other  data  on  waves 
makes  it  possible  to  find  from  (13)  the  dependence  of  the  turbulence  coefficient  on 
the  wave  age  3-  In  view  of  Eqs.  (13)  and  (18),  relations  of  this  kind  with  a = 1 
and  a = 3 in  Eq.  (13)  are  shown  in  Fig.  2 along  with  Neumann’s  empirical  curve. 
Neumann's  curve  was  obtained  by  allowing  for  the  wave  energy  balance.  Thus,  a com- 
parison of  formula  (13)  i^ith  indirect  data  leads  to  the  conclusion  that  a should  be 
close  to  unity.  On  the  other  hand,  the  value  of  exponent  a in  Eq.  (13)  can  be  found 
by  using  expressions  (2,  3,  4)  and  some  hypothesis  on  the  turbulence  scale  & in 
Eq.  (4)  for  the  case  of  wave  motion.  The  general  relations  (7),  (8),  (15),  and  (16) 
which  we  obtained  make  it  possible  to  analyze  some  of  these  hypotheses. 


a > 

r ,*  ‘i 


If, 


f 
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Hypothesis  1.  It  may  be  written  as  follows: 

or  l~  he~  **.  (19) 

T.n  this  case,  using  (2),  (3)  and  (4),  an  expression  corresponding  to  a = 3 in 


Eq.  (13)  is  obtained  for  the  coefficient  of  integral  turbulent  viscosity  V ^ 


As  is  evident  from  Fig.  2,  the  form  of  the  function  $(B)  is  inconsistent  with  con- 
dition (16)  in  this  case.  Thus,  hypothesis  (19)  inadequately  describes  the  turbu-  / 7 76 
lence  conditions  in  a wave,  since  it  leads  to  the  conclusion  that  as  the  scale  of 
wind-driven  waves  increases,  the  turbulence  coefficient  remains  constant. 

Hypothesis  2*  is  related  to  the  assumption  that 

£ - A.  (20) 

It  is  easy  to  see  from  Eq.  (3)  that  this  case  corresponds  to  the  Karman  simi- 
larity hypothesis,  according  to  which 


0u  I 

1 ~ di  I di*  ■ 


(2») 


Using  (2),  (3),  (4),  and  (20),  one  can  obtain  for  V .an  expression  corresponding 

turb 

to  a = 1 in  Eq.  (13).  It  is  evident  from  Fig.  2 that  this  case  satisfies  condition 
(16).  However,  the  case  a = 1,  as  can  be  readily  seen  from  formulas  (7),  (9),  and 
(10),  leads  to  violation  of  condition  (8),  since  when  a = 1,  in  formula  (7) 

ip(6)  = const.  (22) 

The  significant  deficiencies  of  hypothesis  (20)  are  chiefly  due  to  the  following. 
Near  a free  surface,  the  vertical  turbulence  scale  cannot  be  greater  than  the  wave 
height.  Present  concepts  and  data  on  turbulent  fluctuations  at  a free  surface  indi- 
cate their  fairly  high  isotropy.  Thus,  the  only  assumption  would  be  to  consider 
that  near  a free  surface,  the  mixing  length  should  be  proportional  to  the  wave  height, 
not  the  wave  length.  In  addition,  the  free  surface  has  a damping  effect  on  tur- 
bulent fluctuations.  This  effect  is  somewhat  analogous  to  that  of  a solid  wall, 
namely,  with  increasing  distance  from  the  free  surface,  the  turbulence  coefficient 
should  increase  in  some  layer.  As  will  be  seen  below,  this  fact  is  confirmed  by 
certain  indirect  empirical  data.  However,  in  the  case  of  (20)  and  (21),  such  a 
layer  cannot  be  obtained,  since  what  then  occurs  is  only  the  attenuation  of 
turbulence  from  the  physical  liquid  surface  itself,  together  with  the  attenuation  of 
wave  motion.  In  this  connection,  Levich^  proposed  the  following  hypothesis  of  the 
turbulence  scale  in  wave  motion: 


Hypothesis  3.  The  quantity  £ in  Eq.  (4)  is  proportional  to  the  wave  height 
and  increases  with  the  distance  from  the  free  surface  according  to  the  expression 


£ ~ h + z. 


(23) 


Using  (2),  (3),  (4),  and  (23),  one  can  obtain  (see  Ref.  14)  the  following  ex- 
pression for  the  coefficient  V.  . : 

* turb 


Such  a case  was  discussed  In  Refs.  7 and  11.  Hypothesis  (20)  was  first  used 
by  Dobroklonskiy,^ ^ hut  the  results  he  obtained  do  not  pertain  to  the  characteristics 

of  integral  turbulent  viscosity  since  he  did  not  calculate  the  energy  dissipa- 

tion in  trochoidal  waves.  For  this  reason,  the  results  obtained  in  Ref.  12  are  not 
analyzed  in  the  present  paper. 
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hX  . ~ 


where  is  some  constant  coefficient,  and 


;(S)=--(l  + 3*8  + -2^8*). 


We  see  that  hypothesis  (23)  leads  to  results  consistent  with  general  formulas 
(7)  and  (8),  in  contrast  to  hypothesis  (20).  It  is  easy  to  see  from  Eq.  (25)  that 
if  we  neglect  the  first  two  and  last  two  terms  in  Eq.  (25),  we  obtain  the  results 
flowing  correspondingly  from  hypotheses  (19)  and  (20),  whose  deficiencies  we  have 

already  noted.*  On  the  other  hand,  it  follows  that  when  Eq.  (24)  and  (25)  are  used, 

the  dependence  v = $(3)  will  be  represented  in  Fig.  2 by  a curve  occupying  an 

intermediate  position  between  the  curves  a = 1 and  a = 3,  i.e.,  will  satisfy  the 
requirement  formulated  above  for  the  form  of  the  function  <}>(3)  (13),  and  will  be 
similar  to  Neumann's  empirical  curve.  In  addition,  it  follows  from  expressions  (3), 
(4)  and  (23)  that  in  a certain  layer  (z  < X),  the  turbulence  coefficient  will  in- 
crease with  depth.  This  result  may  be  obtained  in  another  way,  which  is  of  interest, 
since  expression  (23)  does  not  flow  directly  from  dimensional  analysis  (in  wave 
motion,  three  parameters  have  the  dimension  of  length:  h,  X,  and  z) . The  transfor- 
mations (9)-(14)  which  we  made  showed  that  the  turbulence  coefficient  may  always 
be  represented  as  a function  of  three  parameters: 

V*  *•  P>-  (26) 

One  can  consider  waves  with  fixed  3 (at  a given  stage  of  development)  and  /' 

assume  that  the  turbulence  coefficient  also  depends  on  the  vertical  coordinate  z. 

Then 

e-  *)•  <27> 

In  view  of  (15),  dimensional  analysis  makes  it  possible  to  write,  by  (27), 


W3  . / g*  \ 
® \W3)  • 


tvrb  g 


It  is  evident  from  Eq.  (28)  that  when  analyzing  processes  of  turbulent  exchange 
in  the  surface  layer  of  the  sea,  which  is  included  in  the  waves,  one  must  consider 
the  dimensionless  parameter  gz/W^  (Froude  number) . This  distinguishes  the  problem 
of  turbulence  in  the  upper  layer  of  a sea  of  homogeneous  density  from  the  analogous 
problem  for  the  ground  layer  of  the  atmosphere. 


In  the  case  z < — , one  can  confine  oneself  to  the  first  terms  of  the  power 
8 / 7 \ 

series  expansion  of  the  function  <!>(  ---  ) and  set  in  Eq.  (28) 

WJ 

£ 

Levich^  assumed  that  the  first  two  terms  in  Eq.  (25)  are  small  in  comparison 
with  the  third.  Actually,  they  are  rather  large,  but,  as  we  have  shown,  strictly 
speaking,  only  one  first  term  in  Eq.  (25)  cannot  be  left  either. 
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® 1 ^ jy*  • 


(29) 


where  A is  some  constant,  which  in  the  general  case  may  be  regarded  as  a function 
of  8.  Comparing  Eq.  (29)  with  the  experiment,  we  will  indicate  the  paper  of 
Pivovarov,15  who,  on  the  basis  of  data  on  propagation  of  diurnal  temperature  waves 
in  the  upper  layer  of  the  sea,  obtained  a linear  dependence  of  the  mixing  coeffi- 
cient on  depth,  and  found  that  A > 0.  Unfortunately,  no  data  are  given  in  this 
paper  on  the  wave  and  wind  regime  for  this  case.  Such  data  are  given  in  Ref.  4, 
but  the  latter  does  not  explain  for  what  level  the  values  of  the  mixing  coefficient 
were  calculated.  Therefore,  it  should  be  emphasized  at  this  point  that  the  rela- 
tions (28)  and  (29),  obtained  from  dimensional  considerations,  together  with  the 
indirect  data  of  Ref.  15  confirm  hypothesis  (23),  according  to  which  the  turbulence 
coefficient  should  increase  in  a certain  surface  layer  with  increasing  distance 
from  the  free  surface  of  the  liquid.  Obviously,  below  this  layer,  because  of  the 
rapid  attenuation  of  the  wave  motion  itself  with  depth,  the  turbulence  should  also  be 
attenuated.  Thus,  the  presence  of  a maximum  of  the  turbulent  mixing  coefficient  at 
some  intermediate  depth  in  the  upper  layer  of  the  sea  is  physically  more  justified. 

In  the  general  case,  the  position  of  this  maximum  should  depend  on  the  wind  velo- 
city and  stage  of  wave  development  8.  These  facts  are  essential  in  solving  certain 
problems  of  turbulent  mixing  in  the  upper  layer  of  the  sea. 

In  conclusion,  the  author  considers  it  his  duty  to  express  his  thanks  to  V.  G. 
Levich  for  his  counsel. 
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